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For the spin models with continuous symmetry on regular lattices and finite range of interac- 
tions the lower critical dimension is d=2. In two dimensions the classical XY-model displays 
Berezinskii-Kosterlitz-Thouless transition associated with unbinding of topological defects 
(vortices and antivortices) . We perform a Monte Carlo study of the classical XY-model on 
Sierpihski pyramids whose fractal dimension is D = log4/log2=2 and the average coordination 
number per site is ss 7. The specific heat does not depend on the system size which indicates 
the absence of long range order. From the dependence of the helicity modulus on the cluster 
size and on boundary conditions we draw a conclusion that in the thermodynamic limit there 
is no Berezinskii-Kosterlitz-Thouless transition at any finite temperature. This conclusion is 
also supported by our results for linear magnetic susceptibility. The lack of finite temperature 
phase transition is presumably caused by the finite order of ramification of Sierpihski pyramid. 
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1. Introduction 

One of the powerful predictions of the renormalization group theory of critical 
phenomena is universality according to which the critical behavior of a system 
is determined by: (1) symmetry group of the Hamiltonian, (2) spatial dimension- 
ality d and (3) whether or not the interactions are short-ranged [1]. The possi- 
bility of phase transitions on systems with nonintegral dimensionality d was first 
considered by Dhar 0] for the classical XY-model and for the Fortuin-Kasteleyn 
cluster model on the truncated tetrahedron lattice with the effective dimensional- 
ity 21og3/log5=1.365. No phase transition at any finite temperature was obtained. 
Subsequently in a series of papers Gefen et al. [3Ha] examined the critical properties 
of the Ising model (discrete Z2 symmetry) on fractal structures, which are scale 
invariant but not translationally invariant, in order to elucidate the relative im- 
portance of multiple topological factors affecting critical phenomena. They found 
that the Ising systems with given fractal dimension D have transition temperature 
T c =0 if the minimum order of ramification R m in > which is the minimum number of 
bonds that needs to be cut in order to isolate an arbitrarily large bounded cluster 
of sites, is finite. In the case of fractals with R m i n = 00 they presented arguments 
that T c >0 for the Ising model. Monceau and Hsiao [?J further studied the Ising 
model on fractals with R m in = 00 using the Monte Carlo method and found weak 
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universality in that the critical exponents depended also on topological features of 
fractal structures. 

For the models with continuous 0(n) symmetry, n >2, on fractal structures 
Gefen et al. [a, 0] used a correspondence between the low-temperature properties 
of such models and pure resistor network connecting the sites of the fractal to argue 
that there is no long-range order at any finite temperature if the fractal dimension 
D <2 even in the case of infinite order of ramification. Subsequently Vallat et al. 
[s[ used the harmonic approximation to the XY- model (0(2) symmetry) on two- 
dimensional Sierpihski gasket (R m i n =3, D = log3/log2=1.585) to show that the 
energy of a vortex excitation is always finite and hence there is no Berezinskii- 
Kosterlitz-Thouless (BKT) transition [it] at any finite temperature as free vortices 
are always present. These conclusions were confirmed in a recent Monte Carlo study 
of the full XY- model on two-dimensional Sierpihski gasket [§] . 

Here we present a Monte Carlo study of the XY-model on three-dimensional 
Sierpihski pyramids with R m i n =A and fractal dimension D = log4/log2=2. The 
model is described by the Hamiltonian 

F = -J^cos(0i-0j), (1) 

where 0< Oi <2tt is the angle/phase variable on site i, denotes the nearest 

neighbors and J >0 is the coupling constant. In the case of translationally invariant 
system in two dimensions this Hamiltonian gives rise to BKT transition and we 
investigate if the same holds in the case of a fractal structure with fractal dimension 
D =2. 

The rest of the paper is organized as follows. In Section 2 we present our algo- 
rithm for generating Sierpihski pyramids and outline the Monte Carlo procedure 
for calculating the thermal averages. Section 3 contains our results and discussion 
and in Section 4 we give a summary. 



2. Numerical procedure 

The procedure which we used to generate three-dimensional Sierpihski pyramids 
(SP) is illustrated in Figure 1, which shows the transition from the zeroth-order SP 
(the tetrahedron of unit side) to the first order pyramid via translations by three 
nonorthogonal basis vectors ei=(l,0,0), e2=(0,l,0) and e3=(0,0,l). The pyramid 
of order m+1 is obtained from the pyramid of order m via translations by vectors 
2 m ei, 2 m e2 and 2 m e3 (m=0,l, ...). It is clear that the number of vertices N m 
of the mth order Sierpihski pyramid can be obtained from the recursion relation 
N m =4N m -i-6 with iVo=4. Thus, in generating the pyramid of order m+1 from the 
pyramid of order m not every point of the mth order pyramid gets translated by all 
three translation vectors 2 m ei, 2 m e2 and 2 m e^. The top of the mth order pyramid, 
(0,0,0), is never translated. The remaining N m -1 points are then all translated by 
2 m ei. Next, the same points except for (2 m ,0,0) are translated by 2 m e2 and finally, 
all points except for (2 m ,0,0) and (0,2 m ,0) are translated by 2 m e3. For pyramids 
of order m <9 we found it most efficient to represent a vertex (i,j,k) by the 
number P = il0 6 +j'10 3 +A:10 and the three translation vectors 2 m ei, 2 m e2 and 
2 m e 3 by the numbers T x (m) =2 m xl0 6 , T 2 (m) =2 m xl0 3 and T 3 (m) =2 m xl0°, 
respectively. Then the result of translating a point represented by number P by a 
vector represented by Tj(m) is described by the number P + Tj(m). From a given 
number P representing a vertex it is easy to get its coordinates k) in the basis 
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{ei, e2, 63}: i = [P/10 6 ], where [• • • ] denotes the integer part, j = [(P— ixl0 6 )/10 3 ] 
and k = P-ixlO 6 - jxlO 3 . 

In the Metropolis Monte Carlo scheme of calculating the statistical averages for a 
model with only the nearest neighbor interactions it is necessary to provide the list 
of nearest neighbors for each site/vertex. We take that all sites that are at a unit 
distance (the size of the edge of the elementary tetrahedron) from a given vertex are 
its nearest neighbors. Thus the number of nearest neighbors of site (i,j,k) varies 
with the order m of Sierpihski pyramid. For example, the vertex (1,1,0) has six 
nearest neighbors in the first order pyramid (Figure 1) but the same site has eight 
nearest neighbors in all higher order pyramids, the additional two neighbors being 
vertices (2,1,0) and (1,2,0). For the Sierpihski pyramids of orders m =4,5,6 we found 
the average number of neighbors per site to be 6.923, 6.981, 6.995 with standard 
deviations 0.903, 0.876, 0.869, respectively. Thus, the average coordination number 
for three-dimensional Sierpihski pyramid is greater than the coordination number 
for three-dimensional simple cubic lattice. This fact should be kept in mind when 
we discuss our numerical results in the next section. In constructing the list of 
nearest neighbors for the sites in a pyramid we found it convenient to group all 
sites according to the values of the sum i + j + k of their coordinates (i,j,k). The 
sites with the same i + j + k belong to the same plane parallel to the basal plane of 
the pyramid defined by points (2 m ,0,0), (0,2 m ,0), (0,0, 2 m ), where m is the order of 
the pyramid (see Figure 1). The nearest neighbors of a given site are located in the 
plane to which it belongs and in the neighboring planes. Throughout this work we 
employed two types of boundary conditions: closed, where the four corners of an 
mth order pyramid were considered to be coupled to each other and open, where 
the four corners are uncoupled to each other. 

The Monte Carlo (MC) simulation of the classical XY-model on Sierpihski pyra- 
mids was based on Metropolis algorithm [l(| • We considered the pyramids of orders 
m =4 (130 sites), m =5 (514 sites) and m =6 (2050 sites). For a pyramid of given 
order the simulation would start at a low temperature with all phases aligned. 
The first 120,000 steps per site (sps) were discarded, followed by seven MC links 
of 120,000 MC sps each. At each temperature the range over which each angle 6% 



was allowed to vary [111] was adjusted to ensure an MC acceptance rate of about 
50%. The errors were calculated by breaking up each link into six blocks of 20,000 
sps, then calculating the average values for each of 42 blocks and finally taking 
the standard deviation a of these 42 average values as an estimate of the error. 
The final configuration of the angles at a given temperature was used as a starting 
configuration for the next higher temperature. 



3. Numerical Results and Discussion 

The heat capacity per site shown in Figure 2 was calculated from the fluctuation 
theorem 

1 (ff 2 ) - {Hf 
C ~ N k^ ' (2j 

where kg is the Boltzmann constant, T is the absolute temperature and (• • • ) 
denotes the MC average. The results did not depend on the type of boundary 
condition (closed or open) within the error bars in analogy to what was found 
for two-dimensional Sierpihski gasket [9j. In the same Figure we also show the 
results for heat capacity obtained for the XY-model on three cubic lattices with 
the periodic boundary conditions. The sizes of cubic lattices were chosen so that 
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they are comparable to the sizes of three Sierpihski pyramids. The peak in the 
specific heat of cubic lattices near ksT / J =2.2 increases with system size fl~2l ] as 
a consequence of diverging correlation length at a continuous (i.e. second order) 
phase transition. On the other hand the specific heat for the Sierpihski pyramids is 
virtually size-independent indicating the absence of long range order for the XY- 
model on three-dimensional Sierpihski pyramid at any finite temperature. Thus, 
although the average coordination number for the Sierpihski pyramid (~ 7) is 
higher than its value for the cubic lattice, the topological properties of this fractal 
structure, in particular a finite order of ramification [5], are responsible for the 
lack of long range order at finite T. The regular lattices have an infinite order 
of ramification as the number of bonds one needs to break in order to isolate an 
arbitrary large bounded cluster of lattice sites is infinite. For fractals with a finite 
order of ramification an arbitrarily large bounded cluster can be cut off from the 
rest of the structure by breaking off only a finite number of bonds and at finite 
temperature thermal fluctuations are sufficient to destroy the long range order. 

The absence of size dependence of the peak in C is found for the XY-model in two 
dimensions 13J, llJ]. In that case the peak results from unbinding of vortex clusters 



131 ] with increasing temperature above the Berezinskii-Kosterlitz-Thouless (BKT) 



transition t emp erature T c where the heat capacity has an unobservable essential 
singularity [151 ]. However the results in Figure 2 for fractal structures of fractal 
dimension D =2 do not necessarily imply BKT transition since size-independent 
peaks in C could result from the average energy per site (E) changing monotonically 
from the values near -3.5J (each site has about 7 neighbors on average) at low 
temperatures to near zero in high temperature paramagnetic phase. 

The main signature of BKT transition is a universal jump in helicity modulus 
7(T c )/r c =2/7r [16[] at critical temperature T c . The helicity modulus "y(T) measures 
the stiffness of the angles with respect to a twist at the boundary of the system. 
At zero temperature, when the angles are all aligned, one expects a finite value of 
7 and at high temperatures, when the system is in disordered paramagnetic phase, 
7 vanishes. For the XY-model on three-dimensional regular lattices j(T) decreases 
continuously with increasing temperature and just below the transition tempera- 
ture T c it obeys a power law 7 (T) oc \T-T C \ V 0. In two dimensions, however, there 
is a discontinuity in 7 stemming from unbinding of the vortex-antivortex pairs at 
BKT transition from quasi-long-range order (order parameter correlation function 
decays algebraically) to disordered phase (order parameter correlation function 
decays exponentially). In numerical simulations on finite systems discontinuity is 
replaced by continuous drop in j(T) which becomes steeper with increasing system 
size (see, for example, [§]). We should point out that since Sierpihski pyramid is 
three-dimensional object topological defects are not only vortices in planes parallel 
to the faces of the pyramid (which are two-dimensional Sierpihski gaskets) but also 
vortex strings. Kohring el al. fl7t ] presented Monte Carlo evidence that a continu- 
ous phase transition for the XY-model in three dimensions is related to unbinding 
of vortex strings. 

We computed the helicity modulus for Sierpihski pyramids following the proce- 
dure of Ebner and Stroud [18] . The idea is to think of Hamiltonian (TjQ) as describing 
a set Josephson coupled superconducting grains in zero magnetic field, where Oi 
is the phase of the superconducting order parameter on grain i. Then if a uni- 
form vector potential A is applied the phase difference Oi — Qj in ([1]) is shifted by 
2-7rA • (rj — rj)/$0) where is the position vector of site i and $0 = hc/2e is the 
flux quantum. The helicity modulus is obtained from the Helmholtz free energy F 
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as 7 = (d 2 F/dA 2 ) A=0 , i.e. 

7 = ( ( dl2 ) A=0 ) " k^f ( ( al) A=0 ) + k^f( 



dH 
~dA 



(3) 

A=0 ' 



We took A to be along one of the edges of elementary tetrahedron (Figure 1). Our 
results for j(T) obtained with closed boundary condition are shown in Figure 3. 
They are completely analogous to what was obtained for two-dimensional Sierpihski 
gaskets @): a rapid downturn in ~f(T) starts near the universal 2/7r-line but the low 
temperature values of 7(T) decrease with increasing system size and the onset 
of the downturn, which is in the vicinity of the putative phase transition, shifts 
to the lower temperatures. For the XY-model on the square lattices, where BKT 
transition does occur, the low-T values of helicity modulus and the onset of its 
downturn do not depend on the system size (see, for example, Q). Our results 
in Figure 3 suggest that in thermodynamic limit N — > oo y(T) vanishes at any 
temperature T >0 implying no BKT transition for the XY-model on Sierpihski 
pyramid of fractal dimension D =2. This conclusion is reinforced by our results 
for 7(T) obtained with the open boundary condition when four corners of an mth 
order pyramid are not coupled to each other, Figure 4. The helicity modulus is 
zero within the error bars which are larger than those obtained with the closed 
boundary condition. The results in Figures 3 and 4 indicate that closed boundary 
condition introduces additional correlations as was the case for two-dimensional 
Sierpihski gaskets 0]. 

Our conclusion about the lack of finite temperature BKT transition is supported 
by results for linear susceptibility 

X N 2 k B T ' { 1 

where M is the magnetization of the system, shown in Figure 5. For finite cubic 
lattices one gets a peak in x near fcgT/J =2.2 whose size and sharpness increase 
with the number of sites as a result of diverging correlation length at the onset of 
long range order. In the case of Sierpihski pyramids the peak in \ a ls° grows with 
increasing system size but it also shifts substantially to lower temperatures. For 
BKT transition Kosterlitz predicted [19] that above T c the susceptibility diverges 
as x ~ exp[(2 — rj)b(T/T c — l) _y ], with 77=0.25, b sal. 5 and z^=0.5, and is infinite 
below T c . Our results suggest that in thermodynamic limit N — > oo there would 
be no divergence in \ at any finite temperature for the classical XY-model on 
Sierpihski pyramid. 



4. Summary 

From our Monte Carlo simulation results we conclude that there is no finite tem- 
perature phase transition for the classical XY-model (0(2) symmetry) on three- 
dimensional Sierpihski pyramid (fractal dimension D=2). Since the heat capacity 
per site does not depend on the system size there can be no long range order at 
any finite temperature. Because the low-temperature helicity modulus decreases 
with increasing system size for closed boundary condition, and is zero within the 
error bars for open boundary condition, it must vanish in thermodynamic limit at 
any finite temperature. This implies no continuous finite temperature phase tran- 
sition associated with unbinding of vortex strings [17( in which case the helicity 
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modulus vanishes at transition temperature T c as a power law j(T) oc \T — T c \ v 
121 ]. Moreover there is no finite temperature Berezinskii-Kosterlitz-Thouless tran- 
sition associated with unbinding of vortices and characterized by discontinuity in 
7(T) at T c . These conclusions are supported by our results for linear magnetic 
susceptibility. The lack of finite-temperature long range order and the vanishing of 
spin stiffness/helicity modulus are the consequence of finite order of ramification 
of Sierpihski pyramid: as an arbitrarily large bounded cluster of sites can be dis- 
connected by cutting only the finite number of bonds thermal fluctuations drive 
helicity modulus to zero and destroy long range order. 
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Figure 1. The first step in creating Sierpinski pyramids. 
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Ar—A N=130 (Sierpinski pyramid) 
▼ — ▼ N=514 (Sierpinski pyramid) 
O — N=2050 (Sierpinski pyramid) 
N=125 (cubic lattice) 
N=512(cubic lattice) 
B—n N=2197 (cubic lattice) 
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Figure 2. Calculated heat capacity per site as a function of temperature. The top three curves are obtained 
for the classical XY-model on cubic lattices while the results for Sierpinski pyramids fall on the same curve 
within the error bars. 
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Figure 3. The helicity modulus 7 as a function of temperature for Sicrpinski pyramids of different sizes. 
The straight line gives the size of Kostcrlitz-Thoulcss discontinuous jump in 7 at various temperatures. 




Figure 4. The hclicity modulus 7 as a function of temperature for Sicrpinski pyramid of order m=5 (514 
sites) obtained with the open boundary condition. 
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^-A N=130 (Sierpinski pyramid) 
N=514 (Sierpinski pyramid) 
3 — N=2050 (Sierpinski pyramid) 
N=125 (cubic lattice) 
N=512 (cubic lattice) 
b—b N=2197 (cubic lattice) 




Figure 5. The linear susceptibility for for Sierpinski pyramids and cubic lattices. Note that the values 
obtained for the cubic lattices have been increased by a factor of 4 in order to fit the scale set by the values 
obtained for the pyramids. 



